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QUANTUM NATURE OF CYCLOTRON HARMONICS IN THERMAL SPECTRA OF NEUTRON STARS 

V. F. SuLEiMANOv^'^, G. G. Pavlov'', and K. Werner^ 
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ABSTRACT 

Some isolated neutron stars show harmonically spaced absorption features in their thermal soft 
X-ray spectra. The interpretation of the features as a cyclotron line and its harmonics has been 
suggested, but the usual explanation of the harmonics as caused by relativistic effects fails because 
the relativistic corrections are extremely small in this case. We suggest that the features correspond 
to the peaks in the energy dependence of the free-free opacity in a quantizing magnetic field, known as 
quantum oscillations. The peaks arise when the transitions to new Landau levels become allowed with 
increasing the photon energy; they are strongly enhanced by the square-root singularities in the phase- 
space density of quantum states in the case when the free (non-quantized) motion is effectively one- 
dimensional. To explore observable properties of these quantum oscillations, we calculate models of 
hydrogen neutron star atmospheres with B ~ lO^^-lO^^ G (i.e., electron cyclotron energy Ec^e ~ 0.1-1 
keV) and Tcs = 1-3 MK. Such conditions are thought to be typical for the so-called central compact 
objects in supernova remnants, such as IE 1207.4-5209 in PKS 1209-51/52. We show that observable 
features at the electron cyclotron harmonics form at moderately large values of the quantization 
parameter, bcs = Ec^e/kTcs — 0.5-20. The equivalent widths of the features can reach ~ 100-200 eV; 
they grow with increasing feoff and are lower for higher harmonics. 

Subject headings: radiative transfer — stars: neutron — stars: magnetic fields — pulsars: individual 
(IE 1207.4-5209, PSR J1210-5226, PSR J1852-h0040, PSR J0821-4300) 



1. INTRODUCTION 

Observations with the X-ray observatories Chandra 
and XMM-JVewton have led to the discovery of absorp- 
tion features in the thermal spectra of several isolated 
(nonaccreting) neutron stars (NSs). Studying such fea- 
tures, one could infer the chemical composition, mag- 
netic field, and gravitational redshift at the NS surface, 
but identification of the features is rather complicated, 
mostly because of the effects of the a priori unknown 
NS magnetic field on radiative transitions in the X-ray 
range. 

First absorption features in the spectrum of an iso- 
lated NS, centered at about 0.7 and 1.4 keV, were discov- 
ered by Sanwal et al. (2002) in Chandra observations of 
IE 1207.4-5209 (hereafter IE 1207), which is the central 
compact object (CCO) of the supernova remnant (SNR) 
PKS 1209-51/52. Similar to other CCOs (Pavlov et al. 
2002, 2004; de Luca 2008; Halpern & Gotthelf 2010), 
IE 1207 does not show the usual pulsar activity, such as 
radio and 7-ray emission or a pulsar wind nebula. This 
NS has the period P = 0.424 s (Zavlin et al. 2000) and 
shows a thermal- like X-ray spectrum (Toff ^ 1-3 MK, 
depending on the model; Zavlin et al. 1998). Further 
observations with XMM-Newton suggested the presence 
of two more features, at 2.1 and 2.8 keV (Bignami et al. 
2003). 

As the energies of the spectral features are harmon- 
ically spaced, it seems natural to interpret them as 
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(gravitationally redshifted) electron cyclotron line and 
its harmonics (e.g., due to the "resonance cyclotron scat- 
tering", by analogy with accreting X-ray pulsars; Big- 
nami et al. 2003), formed in the magnetic field B = 
6 X 10iO(£:;?°j0.7keV)(l + z) G, where E^^{1 + z) = 
Ec,e — heB/rrieC is the electron cyclotron energy, and z is 
the gravitational redshift. This estimate for the magnetic 
field is consistent with the upper limit Bsd < 3.3 x 10^^ G 
obtained by Gotthelf & Halpern (2007) from the upper 
limit on the pulsar's period derivative P. 

However, the cyclotron harmonic interpretation en- 
counters a serious problem because the ratio of the emis- 
sivities (or oscillator strengths) of consecutive harmon- 
ics is known to be proportional to the parameter = 
m&x{kT,Ec.e)/'mec'^ (see, e.g., Pavlov et al. 1980a, b), i.e., 
the emission and absorption in cyclotron harmonics is an 
essentially relativistic phenomenon. As this parameter is 
quite small for IE 1207, ~ 10"'^, observable cyclotron 
harmonics can hardly be expected, even with allowance 
for the fact that the strength of a feature formed in an op- 
tically thick medium depends on the oscillator strength 
nonlinear ly, and it also depends on other factors, such as 
the temperature gradient at the depth where the feature 
is formed. We should note that the conditions in the rel- 
atively cold, low-field atmospheres of IE 1207 and other 
CCOs (where the free-free absorption and emission dom- 
inate over the electron scattering of photons) are quite 
different from those in accreting X-ray pulsars, where 
the parameter is a factor of 30-100 larger, which ex- 
plains the observability of electron cyclotron harmonics 
in those objects (e.g., two harmonics, in addition to the 
fundamental at 26 keV, have been observed in the binary 
pulsar V0332-h53; Pottschmidt et al. 2005) . 

Resonances at cyclotron harmonics (also known as 
quantum oscillations) can, however, arise even in a non- 
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relativistic plasma due to quantization of rotation of a 
charged particle in a magnetic field. In particular, the 
energy of, e.g., a nonrelativistic electron is 



Pz 

2m« 



(1) 



where 71 = 0, 1,2, ... numerates the Landau levels, and 
Pz is the (continuous) electron's momentum along the 
magnetic field. As a result of the quantization, various 
properties of matter show discontinuities when the num- 
ber of the involved Landau levels changes (e.g., when 
a transition to a higher Landau level becomes allowed 
with increasing the energy of the photon that causes this 
transition). Moreover, because the electron motion is 
continuous only in one dimension, the phase space den- 
sity of final electron states in such transitions experiences 
square-root singularities at these discontinuities, which 
gives rise to logarithmically high resonances at harmon- 
ics of the cyclotron energy (see Section 2 for more detail). 

Such quantum resonances in free-free absorption of 
photons in a magnetic field were first noticed by Pavlov 
& Panov (1976; hereafter PP76) who derived the spectral 
opacities and emissivities of a plasma with kT <^ irieC^ 
for different photon polarizations. As the free-free tran- 
sitions are an important source of opacity in NS atmo- 
spheres, one can expect the resonances to manifest them- 
selves as absorption spectral features in thermal spectra 
of NSs. This was qualitatively demonstrated by Pavlov 
& Shibanov (1978) who considered emission from an "at- 
mosphere" in which the source function in the radiative 
transfer equation was approximated by a linear function 
of optical depth. More realistic, self-consistent NS atmo- 
sphere models have been developed later, starting from 
Shibanov et al. (1992), but they used approximate for- 
mulae for the free-free opacity, in which the quantum 
resonances were neglected. 

Here we present first calculations of realistic magne- 
tized NS atmosphere models with the quantum reso- 
nances at electron cyclotron harmonics in the free-free 
opacity taken into consideration and show that harmon- 
ically spaced absorption features are significant in the 
emergent spectra. These calculations demonstrate that 
thermal spectra of NSs with a surface magnetic field 
~ 10^°-10^^ G can exhibit several such features in the 
observable X-ray range, and, in particular, the features 
observed in the spectrum of 1E1207 can be interpreted 
as due to the quantum effects. 

2. QUANTUM OSCILLATIONS OF OPACITIES 

We consider hydrogen NS atmospheres with magnetic 
fields ^ 10^°-10^"'^ G, such that the electron cyclotron 
energy, i?c,e ~ 0.1-1 keV, and its first several harmon- 
ics are in the soft X-ray range, observable with Chandra 
and XMM-Newton. We assume that the atmosphere is 
hot enough, Tcff > 10® K, to neglect the small fraction of 
neutral atoms. The opacity in such a fully ionized atmo- 
sphere is determined by two processes: free-free absorp- 
tion (inverse bremsstrahlung) and scattering of photons 
by electronf0, modified by the magnetic field. We as- 
sume that the atmosphere is cold enough, Tcff < 10^ K, 

* Scattering by protons, as well as the proton contribution into 
the free-free absorption, can be neglected at the relatively low mag- 
netic fields, when Ec p <^ E. 



to neglect the change of photon energy in the scattering 
process (comptonization effects). Moreover, we calcu- 
late the opacities in the so-called "cold plasma approxi- 
mation" (Ginzburg 1970), which implies that relativistic 
effects, such as the cyclotron emission (hence also ab- 
sorption and scattering) at the harmonics of the electron 
cyclotron energy, can be neglected, as well as the thermal 
motion effects, including the Doppler broadening of the 
cyclotron resonance. 

The transfer of high energy radiation in a magnetized 
plasma is described by two coupled transfer equations for 
the ordinary and extraordinary normal modes (Gnedin 
& Pavlov 1974). The normal mode polarizations and 
opacities depend not only on the photon energy E but 
also on the angle Ob between the wave vector and the 
magnetic field; they can be written as (e.g., Pavlov et al. 
1995) 
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where k^^{E) and fcx are the "classical" free- free and 
Thomson opacities at i? = 0, g±i{E) = g±{E) and go — 
g\\ (E) are the Gaunt factors for radiation polarized across 
and along the magnetic field, 'jr = {2/3)E'^{e^ /hmeC^) 
and ja = 7r [kc\{E) / kT] ga{E) are the radiative and col- 
lisional widths, and el;'^{E, 9b) are the projections of the 
normal mode polarization vectors (see Kaminker et al. 
1982, 1983 for a detailed description of the normal mode 
polarizations and opacities in the cold plasma approxi- 
mation) . 

Various forms for the magnetic Gaunt factors g±,\\{E) 
were derived by PP7(2 who calculated the dependence 
of the Gaunt factors on the dimensionless photon energy 
X = E/kT for several values of the "quantization pa- 
rameter" 6e = Ec^e/kT and obtained simple expressions 
for limiting cases. Examples of the energy dependence 
of the Gaunt factors are shown in Figure 1, for 6e = 5, 
0.05, and 0. In our calculations we used Equations (30) 
and (31) from PP76, 
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^ PP76 used Coulomb logarithms, Aj^ y = (•""/ V3) i instead 
of the Gaunt factors. 
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Figure 1 demonstrates that the energy dependences of 
the magnetic Gaunt factors oscillate around the smooth 
non-magnetic curve, with maxima at the cyclotron en- 
ergy and its harmonics. The oscillations are strong in 
g±{E), especially at larger values of be, being almost im- 
perceptible in g^^{E). If the last term in D is omitted 
(i.e., terms oc me/rrip are neglected), the cyclotron reso- 
nances in g±{E) diverge logarithmically (PP76; see the 
dashed curves in the upper panel of Figure 1). The origin 
of the oscillations and divergences can be understood as 
follows. 

For an electron transition n,pz — >■ n',p'^ with absorp- 
tion of a photon with energy E, the conservation of en- 
ergy is described by the delta-function: 
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where P ~ [2me{E — sEc e+pl/^nie)]^^^ , s = n' — n, and 
9(P^) = 9(i? — s_Ec,e+Pz/2me) is the unit step function. 
The equation for the transition rate contains integrations 
over p'^ and Pz and summations over n and n' (or over 
n and s). Because of the step function in Equation (6), 
a new term in the sum over s is added whenever the in- 
creasing E passes through kEc^e— plf^me {k — 1,2,...). 
Each of the added terms contains a singular factor oc 
{E - kEc,e + pl/2me)~^'^ at E ^ kEc^e - pl/2me + 0. 
Further averaging over p^ with the one-dimensional elec- 
tron distribution function turns the square-root singu- 
larities into logarithmic singularities, oc In |£' — fci?c,e|~^- 
Such behavior is seen in g±{E) but not in g\\{E) because 
in (E) the singularity is canceled by the zero of the 
square of the transition matrix element for the longitu- 
dinal polarization [oc (p'^, —p^)^; see Equations (20) and 
(26) in PP76]. 

There exist a number of mechanisms that broaden the 
peaks E = kE^^e and remove the logarithmic singulari- 
ties (PP76). The last term in D (Equation [5]) accounts 
for one of such mechanisms, associated with the proton 
recoil (the broadened peaks of finite heights are shown 
by solid lines in Figure 1). When the magnetic field 
decreases, the distances between the peaks eventually 
become smaller than their widths, the oscillations dis- 
appear, and the magnetic Gaunt factors turn into the 
Gaunt factor at i? = 0. 

We emphasize that the oscillations in (/_lh {E) are solely 
due to quantization of the electron motion across the 
magnetic field, because of which the motion in only 
one dimension (along the magnetic field) remains non- 
quantized and the square-root singularities appear af- 
ter the integration over the final longitudinal momen- 
tum p'^ (see Equation [6]). M B ^ 0, when the 
quantization disappears, the integration over the final 
(three-dimensional) momentum with allowance for the 
energy conservation, J . . . 5[{p''^ — p^)/2me — E] dP'p' = 
AnnieP J ■ ■ ■ dp' [where p' and p are the moduli of the 
final and initial momentum vectors, and P = {p^ + 
2me£')^^^], does not lead to singularities or oscillations. 

3. ATMOSPHERE MODELS 

To calculate the magnetic NS hydrogen atmo- 
sphere models, we use our recently developed code 
(jSuleimanov et al.] |2009() , assuming the magnetic field 
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Figure 1. Dependence of the transverse and longitudinal Gaunt 
factors on relative photon energy for an electron-proton plasma 
with be = Ec^e/kT = 5 and 0.05. The solid curves show the 
Gaunt factors with proton motion and recoil taken into account, 
the dashed curves show the Gaunt factors in the limit of nip — > oo , 
and the dotted curves are the Gaunt factors at B = 0. 



normal to the stellar surface. We assume full ioniza- 
tion and neglect the vacuum polarization by the mag- 
netic field because relatively weak magnetic fields are 
co nsidered. For the opacities , we use the equations 
by Ivan Adelsberg fc Lai I (|2006( 1 , substituting there the 
above-described magnetic Gaunt factortO; these equa- 
tions reduce to Equation (2) when the proton contribu- 
tion is negligible. 

Assuming the surface gravitational acceleration of 1 x 
10^"' cm s^^, we calculated two sets of the atmosphere 
models. In the first set the magnetic field strength was 
fixed {B = 7 X 10^'^ G) and models with three different 
effective temperatures (Tcff = 1, 1.5 and 3 MK) were con- 
structed. A common property of the emergent spectra, 
as shown in Figure [2] (top panel) , is prominent absorp- 
tion features at the cyclotron energy and its harmonics, 
Ek = kEc,e = O.Slfc keV {k = 1,2,...), which origi- 
nate from the quantum peaks in g±{E). The equivalent 
width Wk of the fc-th feature decreases with increasing 
k (similar to the strengths of the quantum peaks in the 
transverse Gaunt factor) as well as with increasing T^s 
(i.e., with decreasing the effective quantization parame- 
ter, bcff = Ec e/kTes ft! 9.4, 6.3, 3.1 for the three models). 
For instance,' Wi = 240, 200 and 180 eV at T^e = 1, 1.5 
and 3 MK, respectively, while W2 = 190, 140 and 55 eV 



' Ivan Adelsberg fc Lai I I I2006I ) use the magnetic Gaunt factors 
which are applicable only at kT <^ -Bc.e and 



by Nagel (1980) 

E < Ec,e - O(kT) (see the Appendix in Kaminker et al. 1983) 
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Figure 2. Emergent spectra {top) and temperature structures 
(bottom) for NS atmospheres with the magnetic field B = 7 X 10^" 
G (soUd curves) and B = (dashed curves) for three effective 
temperatures, T^g = 1, 1.5, and 3 MK. The dotted curves in the 
top panel show the blackbody spectra for the same temperatures. 



at the same temperatures. 

In the top panel of Figure [5] we also show the emer- 
gent spectra for non-magnetic atmosphere models with 
the same effective temperatures. These spectra are close 
to the "continua" of the corresponding magnetic spectra 
at energies higher than the cyclotron energy, but they 
lie below the magnetic spectra at lower energies as the 
radiation emerges from deeper and hotter layers in the 
magnetic case because of lower opacities at i? ^ £^c,e- 

The second set consists of four models with different 
magnetic field strengths {B = 1, 4, 7, and 10 xlO^" G) 
at the same effective temperature (Tcff = 1.5 MK). The 
emergent spectra are shown in the top panel of Figure 
[HI Again, the equivalent widths of the absorption fea- 
tures decrease with increasing harmonic number k and 
decreasing effective quantization parameter {b^e = 0.9, 
3.6, 6.3 and 9.0 for the four models). 

The radiation spectra in the ordinary (O) and extraor- 
dinary (X) polarization modes are shown in the top panel 
of Figure 4, together with the total (X-l-0) spectrum, for 
the model with T^s = 1.5 MK and S = 7 x G. As 
the X-mode opacity is low a.t E Ec e due to the factor 
E'^/{{E ± E^^e? + 7^] « {E/E^^e? <' 1, the X-mode ra- 
diation emerges from deep, hot layers and dominates at 
these energies. But at the cyclotron resonance the opac- 
ity in the X-mode is strongly increased, and the total 
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Figure 3. Top panel: Emergent spectra for magnetic NS atmo- 
spheres with Tefi = 1.5 MK and different magnetic fields (B = 1, 
4, 7 and 10 xlO^'' G), calculated with magnetic and non-magnetic 
Gaunt factors. For clarity, the spectra for Bio = lOi 7 and 1 are 
shifted along the ordinate axis by factors 10+\ lO+O-^ and lO^'^'^. 
Bottom panel: Temperature structures of the models with different 
magnetic fields; 1 (dash-dotted curve), 4 (dashed curve), 7 (solid 
curve), and 10 xlO^" G (dotted curve). 



spectrum is mainly determined by the 0-mode, whose 
opacity is about the same as at i? = 0. As a result, 
the cyclotron absorption feature is very strong in the X- 
mode spectrum but not in the total spectrum because the 
X-mode flux is small (Pavlov & Shibanov 1978). This ex- 
plains the lack of a strong cyclotron feature in the atmo- 
sphere spectra calculated with the non-magnetic Gaunt 
factors (see Figure 3, top panel). Instead of a strong ab- 
sorption line, we see in these spectra a broad, shallow 
depression at E < Ec.e (m the red wing of the cyclotron 
line), where the X-mode starts to dominate in the emer- 
gent spectrum. In contrast, the transverse Gaunt factor 
is present in the opacities of both modes, so that the 
corresponding absorption features are prominent in the 
outgoing spectrum. 

Spectra of emergent specific intensity for the model 
with Toff = 1-5 MK and B = 7 x 10^° G for various 
angles 9 between the line of sight and the surface nor- 
mal are shown in the bottom panel of Figure 21 The 
quantum oscillations in the intensity spectra are seen at 
any angle; the equivalent width of the corresponding ab- 
sorption features depends on 9 nonmonotonously (e.g., 
Wi ~ 220, 160, 320 and 330 eV for 9 = 1°, 30°, 60° and 
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Figure 4. Top panel: Emergent spectrum of the magnetic NS at- 
mosphere with B = 7x 10^" G and T^g = 1.5 MK (solid curve) to- 
gether with emergent spectra in ordinary (O-mode, dashed curve) 
and extraordinary (X-mode, dotted curve) modes. Bottom panel: 
Spectra of the emergent specific intensity for the same model at 
different angles to the surface normal (indicated near the curves). 

89°, respectively). From the comparison of the spectra 
at different 9 we see that the angular distribution of the 
emergent radiation in the absorption features is different 
from that in the continuum. For instance, the specific 
intensities are more peaked towards the surface normal 
at energies outside the absorption features, while the ra- 
diation is generally more isotropic at the center of the 
features. 

The bottom panels of Figures 2 and 3 show the tem- 
perature structures for the six atmosphere models. The 
temperature structure is only slightly affected by the rel- 
atively low magnetic fields. In particular, the magnetic 
field raises the temperature of upper atmospheric layers 
^-t 3 < 6eff ^ 6, apparently because of cyclotron heat- 
ing, while outside of this range the surface temperature 
becomes lower. 

4. DISCUSSION 

We have presented computations of fully ionized hy- 
drogen atmospheres of NSs for magnetic fields B ~ 10^ - 
10^^ G, such that the electron cyclotron energy, Ec^e ^ 
0.1-1 keV, is within the range of energies where the ther- 
mal emission from isolated NSs is usually observed. We 
have shown that the peaks in the energy dependence of 
the free-free opacity, caused by the quantization of the 
electron rotation around the magnetic force lines, lead to 



absorption features (quantum oscillations) at the elec- 
tron cyclotron energy and its harmonics in the atmo- 
sphere spectra. The quantum oscillations are best ob- 
servable at moderately large values of the quantization 
parameter, 0.5 < bcs ^ 20, when the quantization is sig- 
nificant but the features are not too far in the Wien tail 
of the spectrum. The equivalent widths of the absorption 
features reach ~ 100-200 eV in the examples considered; 
they grow with increasing bgff and are lower for higher 
harmonics. 

At least two absorption features, at energies 0.7 and 
1.4 keV, with equivalent widths ^ 100 eV have been 
observed in the X-ray spectrum of the CCO IE 1207 
(Sanwal et al. 2002). Since the thermal and cyclotron 
energies are essentially nonrelativistic for this object 
{kTcs/meC^ < Ec^e/me(? ^ 10"'^), the features cannot 
be explained as due to the commonly known cyclotron 
processes. On the other hand, as 6cff ^ 3-5, the features 
can be naturally interpreted as caused by the quantum 
oscillations. We note that our calculations support the 
reality of the suspecte d features at 2.1 and 2.8 keV in 
the IE 1207 spectrum (IBignami et al. II2003D . To con- 
firm our interpretation and infer the properties of the 
IE 1207's atmosphere, the phase-dependent spectra of 
IE 1207 should be compared with the model spectra ob- 
tained by integrating the specific intensities over the vis- 
ible NS surface at different rotation phases, for various 
orientations of the rotation and magnetic axes. 

In addition to IE 1207, there exist other NSs with sim- 
ilar magnetic fields and temperatures whose emission 
can be described by our atmosphere models. In par- 
ticular, other CCOs, which likely have magnetic fields 
- 10i°-10ii G (Halpern & Gotthelf 2010, and references 
therein), are expected to show quantum oscillations in 
their soft X-ray spectra. However, conditions for their 
observation may be less favorable. For example, for the 
spin-down magnetic field Bgd = 3.1 x 10^° G of the GGO 
J 1852-1-0040 in the Kesteven 79 SNR (Halpern & Got- 
thelf 2010), the spectrum should show absorption fea- 
tures at the cyclotron energy, — 0.36(1 -I- z)~^ keV 
and its harmonics. However, as tlie spectrum below ^ 1 
keV is strongly absorbed by the ISM (TVh « 1.8 x 10^^ 
cm^^), only high harmonics, whose equivalent widths are 
small, might be detectable. Interestingly, the spectrum 
of this CCO, detected by the XMM-Newton EPIC de- 
tectors (see Figure 3 in Halpern & Gotthelf 2010), shows 
a hint of an absorption feature at ~ 1.3 keV, which 
might be third or fourth harmonic of the cyclotron en- 
ergy, E^^ « 0.33 or 0.26 keV, respectiveljQ. To confirm 
the feature and estimate the cyclotron energy and grav- 
itational redshift, a deeper observation of this CCO is 
required. 

For the other pulsating CCO, J0822-4300 in the Pup- 
pis A SNR, the period derivative has not yet been mea- 
sured, and only an upper limit on the magnetic field, 
Bsd < 9.8 X 10^-^ G, has been estimated by Gotthelf 
& Halpern (2009). Fitting the spectrum with a double 
blackbody model, these authors found an emission fea- 
ture around 0.8 keV at some phases. To explain the fact 

^ Halpern &: Gotthelf (2010) do not mention this feature. We 
note that it is seen in the EPIC pn and MOS spectra, which sup- 
ports the reality of this feature despite its low statistical signifi- 
cance in each of the spectra. 
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that the feature is seen in emission, one has to assume an 
energy source in the upper layers of the atmosphere, such 
as the heat released by accreting matter, which is not in- 
cluded in our models. However, it seems that the spectral 
structure at these energies (see Figure 3 in Gotthelf & 
Halpern 2009) can also be interpreted as an absorption 
feature, perhaps a cyclotron line, centered at ~ 0.9 keV, 
which of course would imply a different continuum. To 
examine this interpretation and look for quantum oscilla- 
tions, the CCO's phase-resolved spectra should be fitted 
with our atmosphere models. 

Another class of NSs with possible quantum oscilla- 
tions in their spectra are radio pulsars with magnetic 
fields 10^°-10^^ G. There are 62 pulsars with such fields 
in the ATNF Pulsar Catalogue^, but none of them has 
been observed in X-rays. As these pulsars are old (the 
youngest one, PSR J1810-1820 with B^d = 4.3 x G, 
has the spin-down age r^d — P/2P — 47 Myr), the bulk 
of NS surface is too cold to be seen in X-rays, but their 
polar caps can be heated up to ~ 1-3 MK by relativistic 
particles and 7-rays generated in the pulsar magneto- 
sphere (see, e.g., Zavlin & Pavlov 2004). Observing the 
quantum oscillations in the thermal spectra of polar caps 
is, however, a challenging task because the luminosity of 
this component is low and it is strongly absorbed by the 
ISM for many of the pulsars. 

In conclusion, we should mention that the opacities 
by PP76, used in our atmosphere models, become inac- 
curate close to the centers of the quantum peaks, espe- 
cially at the core of the fundamental cyclotron line. A 
more accurate consideration, which would include cor- 
rections to the Born approximation and the collective 
(high-density) effects (e.g.. Sawyer 2007), and take into 
account the Doppler broadening of the fundamental res- 
onance (Pavlov et al. 1980a), should broaden the spec- 
tral features and decrease the number of observable ones. 
Also, it would be worthwhile to investigate the role of 
the relativistic corrections to the opacities (including the 
resonances at the cyclotron harmonics). Finally, for the 
comparison with the observational data, it would be use- 
ful to construct atmosphere models with magnetic fields 
inclined to the normal to the surface and integrate the 
emission over (a fraction of) the NS surface. We plan to 
consider these problems in future works. 
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